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ABSTRACT 



We present results from our X-ray data analysis of the supernova remnant (SNR) 
G330.2+1.0 and its central compact object (CCO), CXOU J160103.1-513353 (J1601 
hereafter). Using our XMM-Newton and Chandra observations, we find that the X- 
r~| I ray spectrum of J1601 can be described by neutron star atmosphere models (T°° ~ 

2.5-5.5 MK). Assuming the distance of d ^ 5 kpc for J1601 as estimated for SNR 
O ' G330. 2+1.0, a small emission region of ii ~ 0.4-2 km is implied. X-ray pulsations 

■ previously suggested by Chandra are not confirmed by the XMM-Newton data, and 
, are likely not real. However, our timing analysis of the XMM-Newton data is limited 

by poor photon statistics, and thus pulsations with a relatively low amplitude (i.e., an 

intrinsic pulsed-fraction < 40%) cannot be ruled out. Our results indicate that J1601 is 

a CCO similar to that in the Cassiopeia A SNR. X-ray emission from SNR G330. 2+1.0 is 
00 

I dominated by power law continuum (F ~ 2.1-2.5) which primarily originates from thin 

I filaments along the boundary shell. This X-ray spectrum implies synchrotron radiation 

^\ ' from shock-accelerated electrons with an exponential roll-off frequency froUoff ~ 2-3 

X 10^'' Hz. For the measured widths of the X-ray filaments {D ~ 0.3 pc) and the 

■ estimated shock velocity {vg ~ a few x 10'^ km s^^), a downstream magnetic field B 
~ 10-50 iiG is derived. The estimated maximum electron energy -Emax ~ 27-38 TeV 

^ , suggests that G330. 2+1.0 is a candidate TeV 7-ray source. We detect faint thermal 

^ I X-ray emission in G330. 2+1.0. We estimate a low preshock density uq ~ 0.1 cm~'^, 

which suggests a dominant contribution from an inverse Compton mechanism (than the 
proton-proton collision) to the prospective 7-ray emission. Follow-up deep radio. X-ray, 
and 7-ray observations will be essential to reveal the details of the shock parameters 
and the nature of particle accelerations in this SNR. 
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1. INTRODUCTION 

Since n onthermal X-ray syn chrotron emission was discovered in portions of SN 1006's blast 



wave shock (iKoyama et al. 19951 ). several young supernova remnants (SNRs) now show such strong 
particle acceleration sites in the shock front. The detection of nonthermal X-ray emission in SNRs 
thus provides an excellent opportunity to study the generation of high energy cosmic-rays. Based 
on the archival ASCA data, Torii et al. (2006) discovered that the overall X-ray emission from the 
Galactic shell-type radio SNR G330. 2+1.0 shows a featureless spectrum primarily from a power 
law (PL) continuum (F ~ 2.8). G330. 2+1.0 is thus one of the rare members of Galactic SNRs 
in which X-ray emission is dominated by nonthermal continuum from synchrotron radiation of 
shock-accelerated relativistic electrons: there are only three other Galactic SNRs showing such 
characteristics, G347.3-0.5 (a. k. a. RX J1713. 7-3946, Slane et al. 1999), G266.2-1.2 (a. k. 
a. RX J0852. 0-4622 or "Vela Jr.", Slane et al. 2001), and Gl.9+0.3 (Reynolds et al. 2008). 
Torii et al. (2006) noted the general anti-correlation between the X-ray and radio intensities in 
G330. 2+1.0, which may suggest multiple populations and/or acceleration processes to produce 
spatially separated X-ray and radio emission. The ASCA study of G330. 2+1.0 was, however, 
limited by the low photon statistics (~2000 counts for the entire SNR) and poor angular resolution 
(~3' FWHM). 

We observed G330. 2+1.0 with Chandra to perform a detailed imaging and spectral study of 
the SNR, which c ould not be performed with the low angular resolution detectors on board ASCA 



(jPark et al. 20061 ). In our initial work on the Chandra data, we discover ed a candidate n eutron 



star (CXOU J160103. 1-513353, J1601 hereafter) at the center of the SNR (|Park et al. 20061 ). The 
featureless X-ray spectrum of J1601 is well-described by a black-body (BB) model with T = 5.7 MK. 
The high foreground column (A^h ~ 2.5 x 10^^ cm~^) is consistent with that for SNR G330. 2+1.0, 
supporting their spati al association. Assuming the d istance d ~ 5 kpc to G330. 2+1.0 as estimated 



by the HI absorption (jMcClure-Griffiths et al. 200lh . a small area {R ~ 0.4 km) for the BB region 



is implied. The observed X-ray luminosity is Lx ~ 1 x 10^^ ergs s~^ in the 1-10 keV band. No 
counterpart is found in the optical, IR, and radio bands, and a large X-ray-to-optical flux ratio 
(/i-io kcv//v > 9) is estimated. There is no evidence for long-term time variability (in the 1-7 keV 
band) up to ~10 hr time scales in the light curve of J1601. All these characteristics are typical for 
the peculiar manifestation of neutron stars foun d at the center of several SNRs, dubbed "Centeral 



Compact Objects ( CCOs)" (IPavlov et al. 20021 '). A particularly intriguing aspect of J1601 is its 



possible pulsations ( Park et al. 200^ ). Because of the low photon statistics (~600 counts) and the 



long frame-time (3.24 s) of the Chandra data, the detection of the periodicity {P = 7.48 s) was not 
conclusive {^2a significance). 



We have recently performed a follow-up observation of G330. 2+1.0 and J1601 with XMM- 
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Newton to compensate the low photon statistics and poor time resolution of the Chandra data. 
When combined with the high angular resolution Chandra data, the good time resolution and large 
collecting area of XMM-Newton can help reveal the detailed nature of the SNR and the CCO. We 
report here the results from our analysis of G330. 2+1.0 and J1601 using our XMM-Newton and 
Chandra data. Unfortunately, the XMM-Newton data are significantly contaminated by flaring 
background. Nonetheless, using the available data (Chandra + XMM-Newton), we derive some 
fundamental properties of J1601 and G330. 2+1.0. In §[2l we describe the observations and the data 
reduction. X-ray spectral and timing analyses of J1601 are presented in § [3] and § [H respectively. 
We present the spectral analysis of G330. 2+1.0 in §[3 In §[H1 we discuss implications on the results 
from data analysis of the CCO and SNR. Finally, a summary and conclusions are presented in §[71 



2. OBSERVATIONS & DATA REDUCTION 

We observed G330.2+1.0/J1601 with the European Photon Imaging Camera (EPIC) on board 
XMM-Newton Observatory on 2008-03-20 (ObsID 0500300101). The pointing (RA[J2000] = 16^^ 
01™ 3.^4, Dec[J2000] = -51° 33' 53'.'6) is to J1601 which is positioned at the center of the nearly 
circular X-ray shell of SNR G330. 2+1.0. We chose the small-window mode (4f4 x 4^3 field of view 
[FOV] and 6 ms time resolution) for the EPIC pn to search for pulsations of J1601. We chose the 
full-window mode ('^30' diameter FOV) for the EPIC MOS detectors to study the entire SNR. The 
medium filter was used for all detectors. We reduced the data using the Science Analysis System 
(SAS) software package v7.1.0. 

Our XMM-Newton observations of G330.2+1.0/J1601 were significantly contaminated by flar- 
ing particle background. We removed time bins in which the overall count rate is 2a (the pn) or 3a 
(the MOSl and M0S2) above the mean value for time intervals unaffected by flaring background. 
Time intervals including a considerable contamination by the flaring background (^ 50% above the 
average quiescent rate) were eliminated by these time-fllters. After the time filtering, 26, 31, and 
33 ks exposures for the pr0, MOSl, and M0S2, respectively, are available for further data analysis, 
which is ~40-45% of the total exposure. We then reduced the data following the standard screening 
of event pattern (PATTERN < 12 for the MOSl/2 and PATTERN < 4 for the pn) and hot pixels 
(FLAG = 0). (For the timing analysis of J1601, we used a longer exposure while choosing a smaller 
aperture and more strict event pattern criteria as described in §|4l) There are stable components 
of instrumental background in the EPIC detectors. The primary components that could affect this 
work are Al-K {E ~ 1.5 keV) and Si-K (E ~ 1.7 keV) fluorescence lines due to the interactions of 
high energy particles with the structure surrounding the detectors and the detectors themselve^. 
We removed these events from our image analysis by excluding narrow energy bands centered on 
these lines. Our background-subtracted source spectra show little evidence of these lines. Thus, 



^The ~30% deadtime-corrected exposure for the small window mode of the pn is 18.3 ks. 
^XMM-Newton Users' Handbook, § 3.3.7.2. 
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we believe that the impact of contamination from this instrumental background on our EPIC data 
analysis is negligible. 

Because of the severe contamination by the flaring background, photon statistics of the filtered 
XMM-Newton data are significantly lower than originally intended. Thus, in addition to the XMM- 
Newton data, we use the Chandra data (ObsID 6687) for the spectral analysis to improve overall 
photon statistics. The high angular resolution of Chandra data is also essential to measure the 
widths of the thin X-ray filaments of G330. 2+1.0. We performed the Chandra observation of 
G330. 2+1.0 with the I-array of the Advanced CCD Imaging Spectrometer (ACIS) on 2006-05-22 
as part of the Guaranteed Time Observations program. The effective exposure after the data 
screening is ~50 ks, and thus photon statistics for the SNR and CCO in the ACIS data are similar 
to those obtained by the EPIC M0S1+M0S2 data. The details of the Chandra observation and 
data reduction are described by Park et al. (2006). 



3. X-Ray Spectrum of the Central Compact Object 

We extracted the spectrum of J1601 (~530, 290, and 310 counts in the 0.5-10 keV band for the 
pn, MOSl, and M0S2, respectively) from a circular region with a radius of 15". The background 
spectrum was estimated from two nearby source- free regions with a radius of 30". The background 
counts contribute ~15% to the pn spectrum and ~9% to the MOS spectra. The background- 
subtracted, deadtime-corrected count rates (in the 15" radius aperture) are ~0.025 (pn) and ~0.009 
counts s~^ (MOSl/2). The Chandra spectrum of J1601 was extracted from an ~2" circular region. 
The background spectrum was extract ed from the surro unding annular region with the inner and 



outer radii of 4" and 15", respectively (jPark et al. 20061 ). The background-subtracted ACIS count 



rate is ~0.012 counts s~^. The total source counts combining all the XMM-Newton EPIC and 
Chandra ACIS data are ~1700 counts, which is about three times higher than those used in the 
previous work. Each source spectrum was binned to contain a minimum of 20 counts per energy 
bin. 

We simultaneously fit four spectra of J1601 obtained by the XMM-Newton pn, MOSl, M0S2, 
and the Chandra ACIS. Initially we fit the spectrum with a BB model. The best-fit BB temperature 
and the absorbing column (Tbb = ^-^toi MK, Nu = 2.46t[]J^ x 10^^ ^^-2^ ^2^^ ^ 73.4/74, errors 
are at 90% confidence level [C.L.], hereafter) are consistent with those by Park et al. (2006). The 
implied emitting area is small (R ~ 0.44 dr, km, where d^ is the distance to the CCO in units 
of 5 kpc), which is also in agreement with the previous work. The observed flux (/x ~ 1-22 x 
10"^'^ ergs cm~^ s~^ in the 1-10 keV band) is consistent with the previous Chandra results as well. 
Although a PL model may also fit the data, a very steep photon index (P = 5.6^0 4) and a high 
iVn = 5.4+0.6 X 10^^ cm'^ (^^/zy = 93.6/74) are implied (The PL fit is not acceptable with xl 
> 2, when A^h is fixed at 2.5-3.0 x 10^^ cm~^). This PL shape is too soft for typical synchrotron 
emission from the neutron star's magnetosphere, and the fit is statistically worse than that by the 
BB model. Thus, we conclude that X-ray emission of J1601 is consistent with a BB spectrum. 
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Using only XMM-Newton data, we estimate the same flux (/i^iokeV ~ 1-25 x 10 ergs cm ^ 
s~^), which indicates that flux variations in the two years between the Chandra (2006-05-22) and 
XMM-Newton (2008-03-20) observations are neghgible (^5%). 

While the BB model can fit the overall X-ray spectrum of J1601, it may not be physically 
adequate to describe thermal emission from a neutron star. A neutron star is not a perfect BB, and 
likely has an atmosphere whose emission from the outermost H-layer may dominate the observed 
spectrum (e.g., Pavlov & Zavlin 2000). The observed spectrum of the thermal radiation from a 
neutron star's surface is substantially affected by the properties of its atmosphere such as chemical 
composition, magnetic field, gravity, and the energy-dependent opacities (Pavlov et al. 1995 and 
references therein). A typical observational effect may be a higher temperature and a smaller 
emitting area than the "true" values when the spectrum is fitted by a simple BB model (e.g., 
Zavlin et al. 1998; Pavlov et al. 2000). Therefore, taking advantage of the improved photon 
statistics in the XMM-Newton + Chandra data, we fit the observed X-ray spectrum of J1601 with 
a hydrogen neutron star atmosphere model (NSA model in XSPEC, Pavlov et al. 1995; Zavlin et 
al. 1996). 

First, we fit the spectrum of J1601 with a single NSA model. The magnetic fields of CCOs may 
be significantly lower or higher than the "canonical" magnetic field of a pulsar, B = 10^^ G (e.g., 
Pavlov &; Zavlin 2000; Bignami et al. 2003). Based on the lack of a strong absorption feature in the 
observed spectrum of J1601, which would be interpreted as an electron cyclotron line in an NSA 
spectrum, we can only exclude fields in the range of ~ (2 — 8) x 10^^ G. We thus fit our spectrum 
using all three magnetic field values for which the NSA models are available in XSPEC: B = Q 
(applicable for low fields, B < 10^*^ G), 10^^, and 10^^ G. We fix the neutron star mass and radius at 
the canonical values Mns = IAMq and R^s = 10 km, which correspond to the gravitational redshift 
parameter = {I - 2GM^/ R^sC^Yl"^ = 0.766, and vary the effective temperature, distance, and 
A^H- The fits are statistically acceptable for all three magnetic field values. For B = and 10^^ G, 
the redshifted best-fit effective temperatures, = 2.6 and 3.6 MK, respectively, are lower than 
the Tbb, while = 5.7 MK for B = 10^^ q ^bout the same as the BB temperaturei. The 
best-fit Ah values, 3.1, 2.4, and 2.7 xlO^^ cm'^, for B = 0, 10^^, and 10^^ G, respectively, are 
consistent with that for SNR G330.2+1.0 (A'h ~ 2.5-3 x 10^^ cm.-^). The best-fit distances, 24, 
169, and 55 kpc, for B = 0, 10^^, and 10^'^ G, respectively, are unreasonably large for a Galactic 
object. To reconcile them with the distance to the SNR (d ~ 5 kpc; McClure-Griffiths et al. 2001), 
we have to assume that the sizes of the X-ray emitting areas are much smaller than R^g = 10 km 



■^The higher T^ff and the correspondingly larger distance in the B — 10^^ G NSA fit are caused by the fact that the 
high-energy part of the observed spectrum coincides with the low-energy wing of the the gravitationally redshifted 
electron cyclotron feature, centered at ~ 9 keV. Therefore, the high-energy tail of the X-ray spectrum is much softer 
than those in the models with very low or very high magnetic fields. Although we cannot exclude this field value 
based on the observational data available, we note that the models with slightly lower fields {B ^ 2-8 x 10^^ G) 
would not fit the data, while the models with higher fields (e.g., B > 2 x 10^^ G) would yield essentially the same 
results as the B = 10^^ G model. 
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{R ~ 2, 0.3, and 0.9 km for B = 0, 10^^, and 10^^ G, respectively). Thus, although the fits with 
the B = and 10^'^ G NSA models yield lower temperature and larger sizes than the BB fit, the 
observed emission cannot be interpreted as emitted from the entire neutron star surface. 

Based on these results, it is natural to assume that the X-ray emission in J1601 originates 
from small hot spots on the neutron star's surface, such as suggested for other CCOs (e.g., the 
ceo in Cassiopeia A, Pavlov et al. 2000). In this scenario, the observed thermal X-ray emission 
consists of two characteristic components: the hot component from a small region(s) and the cool 
emission from the rest of the stellar surface (e.g., Pavlov et al. 2000). Therefore, we fit the observed 
spectrum of J1601 with two-component NSA models, assuming B = 0, 10^^, or 10^^ G, the same 
for both components. For the soft component, we fix the distance to the CCO and the size of the 
emitting region at o? = 5 kpc and i? = 10 km, respectively, while varying the surface temperature. 
For the hard component, both the distance and surface temperature are varied. The foreground 
column iVn is tied common for both components, and then is fitted. The results are summarized 
in Table [H The X-ray spectrum of J1601 and the best-fit two-component NSA model with a high 
magnetic field {B = 10^^ G) are presented in Fig. [H We note that, although the additional soft 
component is statistically not required, implying only an upper limit on the observed fiux (e.g., 
/i-iokoV < 5 X 10~^^ ergs cm~^ s~^ at 90% C.L.), the two-component model likely represents a 
physically more realistic picture than the one-component model to account for the small hot region 
implied by the single NSA model. Therefore, we hereafter discuss the spectral nature of J1601 
based on the two-component NSA model fits. 



4. Search for Pulsations from the Central Compact Object 

Park et al. (2006) searched for X-ray pulsations from J1601 in the 50 ks Chandra ACIS 
observation (3.24 s time resolution) and reported a marginally significant (at a ~ 2o" level) periodic 
signal {P = 7.48 s, a pulsed fraction frmp ~ 30%). One of the goals of our follow-up XMM-Newton 
observation was to test the significance of the previously reported period candidate, and to search 
for periodicity outside the frequency range accessible with Chandra data (The EPIC pn in the 
small-window mode provides much better 6 ms resolution). However, as discussed in § [21 the 
flaring background hampered our timing analysis of the EPIC pn data. For the timing analysis, we 
performed the data reduction following the methods described in § [21 except that (1) we used 44 
ks of uninterrupted pn data after removing major flares, (2) we extracted photons from a smaller 
circular aperture (8'.'4 in radius) than the standard pointlike source extraction area for the EPIC 
(15" in radius), and (3) we applied a stricter screening by selecting events with PATTERN = 0. 
After these data reduction, we obtained a total of 761 photons (including w 30% background). 
The arrival times of these photons were recalculated to the sola r system barycenter u sing the SAS 



barycen tool. As with the Chandra data, we used the Z"^ test (jBuccheri et al. 19831 ) to search for 
periodicities in the 5 x 10~^ — 80 Hz frequency range. We calculated Zl at 3.5 x 10'' equally 
spaced frequencies, which corresponds to oversampling by a factor of 10 compared to the expected 
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width TgpJ^jj ~ 22 ^uHz of the Z\ peaks, and guarantees that we miss no peaks. The most significant 
peaks we find have Zl = 31.82,31.68, and 30.56 at / = 72.363177(5)Hz, 12.011894(5)Hz, and 
3.006574(5)Hz. However, even for the maximum value of Z\ = 31.82, the corresponding significance 
is low: only 56.5%, for the number of independent trials J\f = /max^span ~ 3.5 x 10^. Therefore, 
most likely these peaks are due to the noise. 

We also calculated for m = 1,2 arou nd the tentative p ulsation frequency (/ = 0.1336185- 
0.1336999 Hz) suggested by the Chandra data (|Park et al. 20061 ^. We find a broad peak, = 10.06 



at / = 0.133661(5)Hz, overlapping with the Chandra peak. However, even a single trial significance 
(assuming that the periodicity found in the Chandra data is real and hence the pulsation frequency 
is known) for the peak in the pn data is only marginal (2.7c7 with the corresponding /p = 24%), 
while for a blind search significance of this peak is negligible. 

Finally, we searched for periodicity in the combined Chandra and XMM-Newton data. Allowing 
for a non-zero period derivative typical for anomalous X-ray pulsars (AXPs), we calculate Zf on 
two-dimensional grid: / = 5 x 10~^ — 0.15 Hz and P = — 3 x 10~^^ s s~^. The maximum value 
of Zf is 29.9 at / = 0.1336615 Hz and P = 5 x 10"^'^ s s~^. Although the frequency corresponding 
to the maximum Zf is consistent with that of the peak found in the Chandra data alone, the 
significance of the peak found in the joint data is very low. Thus, we conclude that the tentative 
7.48 s periodicity reported with Chandra is not confirmed by the XMM-Newton data. 



5. Spectral Analysis of the Supernova Remnant 

XMM-Newton images of SNR G330. 2-1-1.0 are presented in Fig. [2j Since the small-window 
mode was used for the pn, the SNR is detected only on the MOS detectors. The Chandra image 
has revealed that G330. 2-1-1.0 is a s hell-type SNR with enhanced emission in the thin SW and NE 



parts of the shell (jPark et al. 20061 ). Our XMM-Newton images confirm this general morphology. 
We further reveal spectral variations across the SNR: i.e., the E region of the shell is softer than 
other regions (Fig. [3]). Also, there is a faint hard extended feature at ^2' SW from the CCO 
(marked with an arrow in Figs. [2] and [3]). These features were not clearly seen in the Chandra data 
because of their positions in the ACIS-I chip gaps. 

We extracted the spectra from bright portions of the SNR shell in SW and NE (Fig. U]). 
The SW spectrum was extracted from the ~1' x 3' brightest filament in the SW shell, which 
contains ~1700 counts (~25% of them are background) for the M0S1-I-M0S2 data. This SW shell 
contains ~1300 counts (including ~13% background) for the ACIS-I data. The NE spectrum was 
extracted from a circular region (~30" in radius) in the NE parts of the shell. This region contains 
~570 counts (~30% of them are background) and ~540 counts (including ~30% background) for 
the M0S1-FM0S2 and the ACIS-I, respectively. We used the 0.5-10 keV band spectrum for the 
spectral analysis, and the source spectra were binned to contain a minimum of 20 counts per energy 
bin. Because of the non-uniform particle background across the MOS detectors, the background 
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spectrum for the XMM-Newton data was carefully selected for faint extended sources. We chose 
a few background regions close to the source regions while avoiding any detected (by Chandra) 
point sources. We find generally consistent results between the background subtracted XMM- 
Newton spectra and the Chandra data. Thus, we believe that our background estimates for the 
XMM-Newton data are acceptable. 

The X-ray spectra of G330. 2+1.0 extracted from the SW and NE regions are shown in Fig. [5j 
Our Chandra and XMM-Newton data show featureless continuum-dominated spectra for the bright 
SW and NE filaments. For each region, we performed a simultaneous PL model fit for all the 
three spectra obtained by the MOSl, M0S2, and ACIS-I (Fig. [5|). The best-fit parameters are 
presented in Table [2j The high absorbing column for the SNR shell is consistent with that for 
the ceo J1601, supporting the SNR-CCO association. The PL photon indices are typical for 
synchrotron emission from highly accelerated relativistic electrons. Thus, we fit these SW and NE 
shell spectra with the SRCUT model, which describes X-ray synchrotron emission from the shock- 



accelerated electrons that are also responsible for the observed radio counterpart (jReynolds 19981 : 



Revnolds fc Keohane 1999l ;i. We assume th e radio spectral index o, — 0.3 (where Sy oc ") as 



measured from the entire SNR ([Green 200 ll ). The results of the SRCUT model fits are presented 
in Table El 

The spectrally-hard, extended emission feature at ~2' SW of the CCO is faint: we obtain 
~490 counts from this feature (MOSl+MOSl) in which ~60% of the photons are the background. 
There is no evidence for line features, and the X-ray spectrum may be fitted by a PL of F ~ 2 
with A^H fixed at 2.6 x 10^^ cm~^. This overall spectral shape, the filamentary morphology and 
size (about a few pc), and the proximity to J1601 raise an intriguing possibility that this feature 
might be related to the CCO (e.g., the pulsar wind nebula). Alternatively, it could be a part of 
the SNR shell. However, reliable spectral modeling of this faint feature is difficult because of the 
poor photon statistics. Thus, we do not attempt any further analysis or discussion on this feature. 
Follow-up deep X-ray observations are required to reveal the origin of this potentially intriguing 
feature. 

On the other hand, the E parts of the SNR shell are spectrally softer than other regions 
(Figs. [2] and [3]). The E region spectrum is extracted from an ~lf3 x 2' region in the E parts of 
the shell (Figs. [l]and[6|). This region contains ~730 counts (~45% of them are background) for 
the M0S1+M0S2 data. Since the central part of this region falls in the ACIS-I chip gap, we use 
only the XMM-Newton data for the spectral analysis. The best-fit PL photon index for region E 
is significantly steeper (F ~ 4-5, xt ~ 1-3-1.4, depending on the assumed Ah) than those for the 
SW and NE regions. In fact, the PL of F = 2.3 (as an average for SW and NE shell) cannot fit the 
observed spectrum of region E (x^ ~ 1.8-3.2, depending on assumed A^h) because of a soft excess 
emission aX E ^2 keV. This suggests the presence of soft thermal emission i n the E part of the shel l. 



Thus, we fit the region E spectrum with a plane-shock (PSHOCK) model (iBorkowski et al. 200ll ) 



Since t he photon statistics are po or for this faint feature, we fixed the metal abundances at the solar 



values (jAnders &: Grevesse 19891 ). Initially we fit the observed spectrum with a single PSHOCK 
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model, assuming that X-ray emission in region E is entirely thermal in origin. Then, we used a two- 
component model (PSHOCK -|- PL) assuming that there is underlying nonthermal emission as seen 
in the SW and NE filaments. Results from these two model fits are statistically indistinguishable 
{Xu ~ 1-2 for either model fit). The main difference is that the best-fit electron temperature {kT 
= 1.4+°| keV) for the single PSHOCK model appears to be somewhat higher than that for the 
PSHOCK + PL model {kT = 0.7+^ J keV). The best-fit volume emission measure (EM) for the 
two-component model is higher by a factor of ^2 than that for the single PSHOCK model. Since 
the uncertainties of these measurements are large because of poor photon statistics, it is difficult to 
discriminate these modeled parameters. Thus, we assume plausible ranges of the best-fit electron 
temperature and emission measure in the following discussion, based on the these two models: i.e., 
kT ~ 0.7-1.4 keV and EM ~ 0.6-1.4 x 10^6 ^^-3_ ^he results from the PSHOCK + PL model 
fit for the E region are summarized in Table [2l 



6. Discussion 

6.1. Characteristics of the Central Compact Object 

The previous Chandra data analysis of J1601 revealed characteristics of a thermal spectrum, a 
location at the center of SNR G330. 2-1-1.0, a pointlike morphology without any extended nebulosity, 
the absence of cou nterparts at other wavelengths, and a large limit on the X-ray-to-optical flux ratio 



(jPark et al. 20061 ). J1601 also shows no evidence for long-term variability, which is confirmed by 
our new XMM-Newton data showing a constant X-ray flux (/i_iokeV ~ 1-2 x 10~^^ ergs cm~^ s^^) 
over the ~2 yr period. Thus, J1601 is most likely the CCO associated with SNR G330. 2-1-1.0. Park 
et al. (2006) noted that the BB temperature is higher than the surface temperature expected from 
the standard cooling of a young neutron star, and that the estimated emitting area is too small to 
be a neutron star. It was also noted by Park et al. (2006) that the suggested candidate pulsations 
with a long period, if confirmed, would have been typical for an AXP. 

Our results from the XMM-Newton and Chandra data analysis indicate that the hot component 
emission (T^ ~ 2.5-5.5 MK, depending on the assumed B) must originate from a small region of 
i?h ~ 0.4-2 km. The estimated size of the hot region varies depending on the assumed values 
of the magnetic field and the distance to the CCO. Nonetheless, within the ranges of parameters 
that we consider in this work {B = 0, lO^'^, or 10^^ G, and d = 5-10 kpc), the size of the hot 
region is significantly smaller than the canonical size of the neutron star; i.e., the largest area could 
be i?h ~ 4 km, where B = and d = 10 kpc. A small hot region(s) has been suggested in other 
CCOs, probably indicating X-ray emission from a locally-heated region such as the hot polar cap 
(e.g., Pavlov et al. 2000). On the other hand, the estimated surface temperature of the neutron 
star is significantly lower (Tg°° < 1.5 MK) than that of the hot region. According to the standard 
cooling curves of a neuron star (e.g., Tsuruta 1998; Yakovlev &: Pethick 2004), this temperature 
limit corresponds to a lower limit of several 10^-10^ yr for the neutron star's age. This neutron 
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star age is in plausible agreement with the estimated age of SNR G330. 2+1.0 (see §§ 16.21 and I 
and Torii et al. 2006). The overall characteristics such as the low T^, the high T^, and the small 
i?h are c onsistent with thos e found in the prot otype CCOs in Galactic SNRs such as Cas A and 



Vela Jr. ( Pavlov et al. 2000l : IPavlov et al. 2001 



Since the X-ray flux from the small hot region contributes a significant fraction of the observed 
flux (> 50% of the total flux in the 1-10 keV band), the observed X-ray emission from J1601 may 
be expected to pulsate. However, our XMM-Newton data do not show any conclusive evidence for 
pulsations, indicating that the previously suggested pulsations are unlikely real. We note that the 
low photon statistics in the EPIC-pn data are not sufficient to detect pulsations with an intrinsic 
pulsed- fraction /p ^40%. With the combined Chandra and XMM-Newton data, the detection of 
pulsations with /p ^25% is not feasible. Thus, the presence of an X-ray pulsar for J1601 is not 
ruled out by the current data. The neutron star's magnetic field, which would provide critical 
information on the nature of the object, remains unknown. Deep X-ray observations of J1601 are 
required to make conclusive remarks on the nature of J1601 such as pulsations, magnetic field, age, 
and the origin of its X-ray emission. 



6.2. Nonthermal X-Ray Emission of the Supernova Remnant 

Our joint spectral analysis of the XMM-Newton and Chandra data of G330. 2+1.0 shows that 
X-ray emission from the bright filaments of the SNR shell is dominated by a PL continuum. We 
find that this PL spectrum prevail s for the most parts of the SNR, which was also suggested by 



a previous study (jTorii et al. 20061 ). The best-fit PL model for the bright SW and NE regions of 
the shell indicates photon indices of P ~ 2.1-2.5 which are typical for synchrotron emission from 
shock-accelerated relativistic electrons. Although thermal plasma models may also fit the observed 
spectra, the estimated electron temperatures are high (kT ~ 4-5 keV), and low metal abundances 
0.1 solar, Anders & Grevesse 1989) are required. While a thermal origin of X-ray emission 
from the SNR shell may not be completely ruled out by the current data, the estimated plasma 
temperature and abundances appear to be unusual for SNRs. Thus, except for region E 16. 3p . we 
discuss this SNR based on the nonthermal interpretations of X-ray emission. 

According to our SRCUT model fits of the SW and NE filaments, the best-f it exponential roll- 



off frequency, z^roUoff ~ 1.6-3.3 x 10^^ Hz, is relatively high amo ng Galactic SNRs (jReynolds &: Keohane 19991 ). 



while being simil ar to those for SN 10 06 (jBamba et al. 20031 ) and the bright TeV 7-ray emitting 



SNR G347.3-0.5 ([Lazendic et al. 20041 ). If the particle (electron) acceleration is limited by syn- 
chrotron losses, the cutoff frequency corresponding to the maximum electron energy -Emax is fm(loss) 
oc B E'^^^{Ioss). Since £'inax(loss) oc B~2, fm(loss) is independent of B, and depends only on the 
shock velocity: e.g., assuming a strong shock of the compression ratio of > 4 and the shock normal 
perpendicular to B, the cutoff frequency is fm(loss) ^ 3 x 10^^ r/ f| Hz, where V3 is the shock 
velocity in units of 10^ km s~^, and the ratio of the electron scattering mean free path to the 
gyroradius rj > 1 (e.g., Reynolds 1998; Lazendic et al. 2004). As discussed below and in § 16.31 
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the shock velocity appears to be roughly Vg ~ 4000 km s~^ for G330. 2+1.0, and thus we estimate 
fni(loss) ^ 5 X 10^'' Hz. Unless the shock velocity is much higher and/or the particle acceleration is 
inefficient {r] » 1), the estimated z/m(loss) is comparable with the observed froUoff, suggesting that 
the particle acceleration of electrons in G330. 2+1.0 is likely limited by synchrotron losses rather 
than the age of the SNR. The peak frequency of a synchrotron emitting electron is fp = 1.8 x 10^* 
E^B Hz, where B is the postshock magnetic field perpendicular to the shock normal, and Eq is the 
electron energy. For Up --^ 7 x 10^^ Hz (or ~3 keV) representing the typical X-ray photons based 
on the observed spectrum of the nonthermal filaments in G330. 2+1.0, the corresponding electron 
energy is -Eg = 0.62 B~2 ergs. The characteristic synchrotron loss time scale for such electrons can 
then be estimated to be rioss = 630 E~^ B~^ s = 1017 B~2 s. 

We estimate rioss by measuring the widths of the bright nonthermal filaments of G330. 2+1.0 
using Chandra images (Fig. [7]). We construct projected intensity profiles across the bright SW 
filaments by averaging the photon counts (in 4" pixel bins) over the 40" segments along the filaments. 
We fit these 1-D intensity profiles with a Gaussian to estimate the widths of the filaments. We note 
that high resolution Chandra images of bright X-ray synchrotron filaments in young SNRs show 
typical substructures of a broad exponential downstream region and a much steeper flux decay in 
the upstream (e.g., Bamba et al. 2003). G330. 2+1.0 is more distant than other young SNRs (that 
show bright nonthermal filaments), and the X-ray shell is relatively faint, which does not allow us 
to resolve such a substructure. Since the downstream region is observed to dominate the width of 
the filaments, we assume a negligible contribution from the upstream emission in the widths of the 
filaments to measure the downstream widths of the filaments with a simple Gaussian model. The 
measured widths are ~12"-16" (FWHM) which correspond to physical sizes D ~ 0.3-0.4 pc. 
Because of the far distance and faint surface brightness of G330. 2+1.0, our width measurements 
could be an overestimate from superpositions of thinner filaments. Nonetheless, the estimated 
widths are comparable with an average value for the individual filaments in SN 1006 (~0.2 pc, 
Bamba et al. 2003). Therefore, we take our measurements as a first-order estimate, and certainly 
as an upper limit. 

The advection distance of the downstream electrons from the shock is -Dad = Vs rioss r~^, 
where r is the compression ratio in the shock. Since the direct measurements of the shock velocity 
of G330. 2+1.0 are not available, we consider some plausible estimates for the shock velocity based 
on several independent approaches. Assuming an electron-ion temperature equipartition in the 
postshock region, the detected thermal emission of G330. 2+1.0 (region E) implies Vs ~ 1000 km 
s~^ (^ 16. 3p . This value may be considered as a lower limit for Vg, because the assumed temperature 
equilibration between electrons and ions may have not been established in relatively young SNRs 



with Vg ^ several 10 km s (jGhavamian et al. 20071 ). G330. 2+1.0 shows similar characteristics 



(e.g., the SNR age, froUoflf; and the physical width of the nonthermal filaments etc.) to those of 
G347.3-0.5 and SN 1006 in which the shock velocities are high {vs ~ 3000-4000 km s"\ e.g., 
Parizot et al. 2006 and references therein). The ambient density for G330. 2+1.0 (no ~ 0.1 cm~^, 
§ 16. 3p is not unusually high compared with other SNRs (e.g., Bamba et al. 2003). Thus, the 
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actual shock velocity of G330. 2+1.0 is likely higher than Vg ~ 1000 km s~^, perhaps close to Vg 
3000-4000 km s~^. In fact, models predict high shock velocities of Vg ^ 2000 km s~^ for an efficient 
particle acceleration (e.g., Ellison et al. 2000;2004). Although a small sample is used, an empirical 
relationship between t'roUoff and the physical width of the nonthermal filaments D is der ived to 



be froUoff D~'^ = 2.6 x 10^''' r^^^ for several young historical SNRs (jBamba et al. 20051 ). This 
empirical relation suggests an SNR age tsnr ~ 1000-1200 yr for G330. 2+1.0 for the measured 
z^roUoff '~ 2-3 X 10^"^ Hz. The inferred young age and the low ambient density suggest that the 
SNR may be in a free-expansion or an adiabatic phase, or could be in transition between the two. 
Assuming an adiabatic phase, the suggested SNR ages imply Vg ~ 2300-2800 km s~^ for the SNR 
radius of i? = 7.3 pc (see ^ l6.3l for the SNR radius). For a free-expansion phase, Vg ~ 5800-7000 km 
s~^ is implied. These high velocities are consistent with those estimated for young SNRs showing 
an efficient particle acceleration (e.g. Parizot et al. 2006 and references therein). Thus, as a rough 
estimate by averaging several values discussed above, we for simplicity adopt a shock velocity Vg ~ 
4000 km s~^ for G330. 2+1.0. This shock velocity is admittedly not a measurement and thus only 
a crude first-order estimate. (We would allow a factor of ~2 uncertainty in this velocity estimate, 
and within this range, our conclusions as discussed below are not affected.) 

Assuming r ~ 5-8 for an efficient particle acceleration (e.g., Ellison et al. 2007), we esti- 
mate Tioss ~ 350-600 yr for the measured D ~ Dad ~ 0.3 pc, and thus B ~ 14-20 ^G. The 



maximum electron energy can be estimated by -Emax = 2.5 x 10 ^ ^roUoff ^ ^ TeV = 100-144 



1 1 



■B^ j TeV, where B/ .g is the postshock magnetic field in units of /xG (jRevnolds &: Keohane 19991 : 



Lazendic et al. 20041 ). Thus, -Bmax ~ 22-38 TeV (depending on measured I'roiioff) is derived. In 
addition, if we consider a geometrical projection effect in measuring the widths of the nonthermal 
filaments (e.g., the observed width is ~ 4.6 x the actual width assuming a spherical shock with an 
exponential emission profile. Ballet 2006), the estimated B can be a few times higher (~50 /iG). 

The estimated -Emax for G330. 2+1.0 suggests that this SNR is a candidate 7-ray source. For 
instance, the 7-ray emission by the inverse Compton (IC) scattering off interstellar photons can be 
estimated by E-y ~ 5.1 x 10~^^ -E* eV, where Ej is the ave rage final energy of the up-scattered 



photons, and E^ is the typical energy for the seed photons (jTatischeff 20081 ). Using E'^ ~ 7 x 
10"'^ eV for the cosmic microwave background (CMB) and E^ = -Emax ~ 30 TeV, we estimate 
E^ ^ 3 TeV. However, G 330. 2+1.0 is not identified in the H.E.S.S. Galactic plane survey catalog 



(jAharonian et al. 2006al ). It is probably because G330. 2+1.0 is more distant and thus apparently 
fainter than other TeV-bright SNRs (e.g., G347.3-0.5 and G266.2-1.2). The IC to synchrotron fiux 
ratio /ic//syn = SnUri^d/B^ ~ 10 B~q ~ 0.004-0.1 (where the energy density of the seed CMB 
photons C/rad ^ 0-25 eV cm~^) for the plausible range of -B ~ 10-50 //G in G330. 2+1.0. These 
/ic//syn are in fact similar to the observed fTev/fx for SNRs G347.3-0.5 and G266.2-1.2 (e.g., 
Matsumoto et al. 2007 an d references there in). Then, the overall X-ray fiux of /syn ~ 10^^^ ergs 



cm-2 s-^ for G330.2+1.0 (ITorn et al. 20061 ) implies fic ~ lO^^^-lO'^^ ergs cm'^ s'^. The sky 
position of G330. 2+1.0 was at the edge of the H.E.S.S. survey, in which the exposure was short (<5 
hr). Considering the small angular size (~10') of G330. 2+1.0, which is close to the point spread 
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function of the H.E.S.S. (several arcminutes), and the short exposure in the survey, the estimated 
IC flux is likely close to or below the H.E.S.S. detection limit of / ~ 10~^^ ergs cm~^ s~^ at E ^ 
1 TeV (e.g., Aharonian et al. 2005). Thus, if the 7-ray emission from G330. 2+1.0 is dominated by 
the IC process of the same electrons to produce X-ray synchrotron emission, the non-detection of 
G330. 2+1.0 with the current H.E.S.S. survey data is not surprising. A deep search of 7-ray emission 
for G330. 2+1.0 using ground-based TeV telescopes and Fermi (formerly GLAST) is warranted. 

It is notable th at nonthermal X- ray emission in G330. 2+1.0 is generally anti-correlated with 



the radio emission (jTorii et al. 20061 ). Our high resolution Chandra and XMM-Newton images 
reveal that there actually exist radio counterparts for the bright X-ray filaments in SW and NE, 
but the radio emission is faint (Fig. [3]). The brightest radio emission is in the E parts of the SNR, 
where X-ray emission is faint and spectrally soft (Fig. [3]). Thus, the bright radio emission likely 
traces high density regions where soft (thermal) X-ray emission is enhanced. Based on our SRCUT 
model fits. X-ray emission in SW and NE filaments implies the 1 GHz radio flux of 0.7-1.5 x 10~^ 
Jy, while the MOST 843 MHz image of the SNR suggests ~0.1 Jy for these regions (assuming that 
the total 1 GHz flux for the entire SNR is 5 Jy, Green 2001). Although our radio flux estimates are 
crude and should be considered only as an order-of-magnitude approximation based on a simple 
"normalization" of the total image intensity to the area corresponding to the X-ray-bright SW and 
NE filaments, the discrepancy is substantial by three orders of magnitudes, and should thus be real. 
We do not have an immediate answer as to what causes the large difference between the modeled 
and observed radio fluxes corresponding to the X-ray bright filaments. One speculation is that the 
radio spectral index might not be uniform across the SNR. While the overall radio spectrum is 
fitted by a = 0.3, the faint radio filaments corresponding to the bright X-ray shell might have a 
steeper spectrum. For instance, if we assume a plausible range of the observed radio flux ~0. 01-0.1 
Jy for the SW region and vary the radio spectral index in our SRCUT model fit, we obtain a 
best-fit a ~ 0.53-0.66 (x^ = 1-2). These radio spectral indices are not unusual for shell-type SNRs 
( Green 200 ll ). The best-fit roll-off frequencies are high, but are poorly constrained (froUoff = 13^9^ 

+25 
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X 10^^ Hz when the 1 GHz radio flux of 0.1 Jy is assumed, and froUoff = 8^1^ x 10^'' Hz for the 



radio flux of 0.01 Jy). Although the high roll-off frequency, froUoff ~ 10^^ Hz, implies somewhat 
higher estimates for the shock velocity and the maximum electron energy, these changes do not 
make a significant effect on our conclusions presented here. High resolution radio observations with 
a deep exposure would be essential to study the detailed relationship between the X-ray and the 
radio emission in this SNR. 



6.3. Thermal X-Ray Emission of the Supernova Remnant 

In the E region, soft thermal emission is a significant component in the observed X-ray spec- 
trum. The best-fit electron temperature is kT ^ 0.7-1.4 keV, depending on models (§ [5|). The 
best-fit ionization timescale appears to be high {n^t ^ 10^^ cm~^ s) suggesting that the plasma 
could be in collisional ionization equilibrium, but the n^t parameter is not well-constrained because 
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of the low photon statistics. Detecting thermal emission in SNRs in which nonthermal emission 
dominates is critical to reveal the environmental conditions (e.g., ambient density) and the su- 
pernova energetics that should have affected the SNR evolution and the particle acceleration. In 
fact, 0330.2+1.0 is the only example to reveal thermal X-ray emission among the four Galactic 
SNRs which have been known to be dominated entirely by nonthermal X-rays (see §[I]). Therefore, 
although it is difficult to perform a thorough spectral analysis of thermal emission and to draw 
firm conclusions on the nature of the SNR because of the poor photon statistics for the faint ther- 
mal component, we present a brief discussion on some fundamental SNR parameters based on our 
spectral analysis of region E. 

Based on the best-fit volume emission measure (EM = UenjiV, where ne, n-n, and V are the 

postshock electron, proton densities, and the X-ray emitting volume, respectively), we estimate Ue 
1 -i 

~ 0.4-0.5 f~^d^ ^ cm~^ (where / is the X-ray emitting volume filling factor). These postshock 

1 -i 

electron densities correspond to the preshock hydrogen density no ~ 0.1 f~^dr^ ^ cm~'^. In these 
estimates, we assume Ue = 1.2 nn for the mean charge state with normal composition, and nn 
= 4no for a strong shock. We use the emission volume 1/ ~ 4 x 10^^ cm~^ assuming that the 
path-length through region E is comparable to the physical size corresponding to the angular size 
(~2') of region E at d = 5 kpc. Assuming an ion-electron temperature equilibration, the measured 
electron temperature implies a shock velocity of Vg ~ 800 {kT = 0.7 keV) - 1100 {kT = 1.4 keV) 
km s~^. However, equipartition of the electron-ion temperatures may not have been reached, and 
thus the actual shock velocity could be higher than Vg ~ 1000 km s~^, probably by a factor of a 
few (§ 16. 2p . We estimate the SNR radius of i? ~ 5' (the half of the angular distance between the 
the bright SW and NE filaments), which corresponds to the physical distance of ~7.3 pc. Then, 
assuming an adiabatic phase for the SNR, we apply the Sedov solution to derive the SNR age tsnr 
~ 1100 d^ yr (e.g., for Vg ~ 2500 km s~^, § 16. 2p . For a free-expansion phase, the SNR age is also 
derived to be tsnr ~ 1100 yr (e.g., for Vg ~ 6500 km s~^, § 16. 2p . Using a Sedov solution, the 

5 

explosion energy is estimated to be Eq ~ 2-9 x 10^° d^ ergs for rgN^ ~ 1000-2000 yr. 



Summary and Conclusions 



Based on the ASCA data, the overall X-ray emission fr om SNR G330. 2+ 1.0 was suggested to 
be continuum-dominated with no evidence for line features (jTorii et al. 20061 ). The high resolution 
Chandra images subsequently revealed that X-ray emission from this SNR originates primari ly from 
the thin shell with enhanced filaments in the SW and NE parts of the shell (jPark et al. 20061 ) . Park 
et al. (2006) have also discovered the CCO J1601 at the center of the SNR. We performed follow- 
up observations of G330. 2+1.0 with XMM-Newton to investigate the nature of the CCO and the 
SNR. Although our spectral and temporal analyses of J1601 and G330. 2+1.0 are limited by poor 
photon statistics of the XMM-Newton data caused by significant contamination from fiaring particle 
background, we find several important characteristics of these objects utilizing the XMM-Newton 
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and Chandra data. 

The X-ray spectrum of J1601 can be described by two-component neutron star atmosphere 
models. X-ray emission primarily originates from a small hot region {R ~ 0.4-2 km, T ~ 2.5-5.5 
MK). The rest of the neutron star's surface is cooler (R ~ 10 km, T < 1.5 MK), suggesting an 
■^Iq3-4: yj. neutron star based on the standard cooling models. The neutron star atmosphere 
models do not provide useful constraints on the magnetic field of J1601 with the current data. The 
previously suggested pulsations {P ~ 7.48 s) are not confirmed by the XMM-Newton data. These 
characteristics are similar to those found for CCOs in other Galactic SNRs such as Cas A and 
Vela Jr. The spectrally hard, faint nebulosity at ~2' SW from the CCO could be the associated 
PWN, but its true nature is uncertain with the current data because of the poor photon statistics. 
Follow-up deep X-ray observations are required to reveal the detailed nature of J1601. 

Assuming that X-ray emission in the shell of G330. 2-1-1.0 is synchrotron radiation from the 
shock accelerated electrons, the roll-off frequency of froUoff ~ 1.6-3.3 x 10^ Hz is estimated. 
It is difficult to measure the shock velocity with the currently available data. Based on several 
independent approaches, we make a rough estimate of the shock velocity Vg 4000 km (with a 
factor of ~2 uncertainty). Based on this shock velocity and the measured roll-off frequency, we find 
that the particle (electron) acceleration in G330. 2-1-1.0 is likely limited by synchrotron losses rather 
than the SNR age. Using the Chandra images, we measure the widths of the bright nonthermal 
X-ray filaments {D ~ 0.3-0.4 pc). Using these widths and the shock velocity, we estimate the 
synchrotron loss time of rioss ~ 350-600 yr and the magnetic field of i? ~ 10-50 fiG. The maximum 
electron energy is derived to be Ey^ax ~ 22-38 TeV. These electron energies suggest that G330. 2-1-1.0 
is a candidate 7-ray source (up to ~TeV) by the IC scattering of the CMB photons. The non- 
detection of G330. 2-1-1.0 in the current H.E.S.S. survey with a short exposure is perhaps expected, 
because G330. 2+1.0 is more distant and likely a fainter 7-ray source than the bright TeV SNRs 
like G347.3-0.5 and Vela Jr. 

G330. 2-1-1.0 is particularly intriguing because this is the only SNR in which we detect a thermal 
component among the four Galactic SNRs known to be dominated by nonthermal X-ray emission. 
Although the uncertainties are large due to the poor photon statistics, the estimated density (no ~ 
0.1 cm~^) is low, suggesting that 7-ray emission, if it exists, would be dominated by the IC process. 
The detection of 7-ray emission as well as thermal X-ray emission with high photon statistics 
from G330. 2-1-1.0 will be essential to test and constrain models for 7-ray production from shock- 
accelerated particles. Follow-up deep observations with X-ray detectors on board XMM-Newton 
and Suzaku are necessary for a thorough study of thermal X-ray emission. Deep 7-ray observations 
using Fermi and the ground-based TeV telescopes will be critical to reveal the nature of nonthermal 
radiation produced by shock accelerated particles. High resolution radio and X-ray observations of 
G330. 2-1-1.0 with a deep exposure are essential to reveal the origin of the apparent inconsistency 
between the radio and nonthermal X-ray emission, such as the radio spectral index variation across 
the SNR. 
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Tablc 1. Best-Fit Parameters of J1601 from two-component neutron star atmosphere models. 



B 

(10^2 Q) 


Nh 
(10^2 cm-2) 


(10« K) 


(10^ K) 


iih'' 
(ds km) 


/l-lO kev'^ 

(10"^^ ergs cm-2 s"^) 


Li-io keV 

(10^3 di ergs s-i) 





1 

10 


q 9Q+0.53 


< 1.4 

< 1.5 

< 1.5 


^■■^-0.2 

5 5+"-'' 

^•^-0.3 


Q 4+0.1 

*-'-'^-0.3 


1.23±0.06 
1.21±0.06 
1.21±0.06 


1.5 
1.9 
1.8 


80.5/78 
79.9/78 
79.7/78 



Note. — Errors are at 90% confidence. Mns and i?ns are fixed at 1.4 M© and 10 km, respectively. 

^Tg"" and are effective temperatures of the cool neutron star surface and the hot small region, respectively, as 
measured by a distant observer where T°° = QrT and Qt = (l-2GMns/JinsC^)2 . 

''The radius of the hot region scaled by d = 5 kpc. 

'^Tho observed flux in the 1-10 keV band. Assuming a Poisson distribution of the observed photon statistics, 2cr 
statistical errors are quoted. 



Table 2. Best-Fit Power Law Model Parameters for SNR G330.2-M.0. 



Region 


(10^2 cm-2) 


r 


kT 
(keV) 


riot 

(10" cm-3 s) 


EM 
(105« cm-3) 




SW 
NE 
E 


2 fin+"-'"' 

^■OU_().34 
o 04+0.65 

2 4t^+0'^2 

^•*O_0.57 


9 1 q+0.24 
^•^''-0.22 
9 r9+0.40 

2.3 


0.70tJ;l 


> 5 




163.9/137 
51.9/51 
35.3/30 



Note. — Errors are at 90% confidence. For region E, parameters firom the best-fit two compo- 
nent model (plane-shock -|- power law, where F = 2.3 is fixed) are presented. 



Table 3. Best-Fit SRCUT Model Parameters for SNR G330.2-M.0. 



Region 


Nh 
(10^2 cm-2) 


I^rolloff 

(10^^ Hz) 


1 GHz Flux 
(mJy) 


x'/y 


SW 
NE 


9 41+0.17 

9 rQ+0.32 
^.03_Q 28 


3.3it;^ 
i.6it;? 


0.15±0.01 
0.07±0.01 


163.7/137 
51.9/51 



Note. — Errors are at 90% confidence. The radio spectral index 
is fixed at a = 0.3 (where S^, oc v'"). 
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Fig. 1. — The X-ray spectrum of J1601 as observed by XMM-Newton and Chandra. The best-fit 
two component NSA model with B = 10^'^ G is overlaid. The solid lines are the best-fit model. 
The dotted- and dashed-lines are the small hot region and the large cool surface components, 
respectively. The lower panel is the residuals from the best-fit model. 
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Fig. 2.— XMM-Newton images (M0S1+M0S2) of G330.2+1.0: (a) the soft band (1-2 keV), and 
(b) the hard band (2-8 keV). In (a), the 1.4-1.6 and 1.7-1.8 keV bands are excluded to remove the 
bright instrumental lines at E ^ 1.5 (Al K) and 1.74 (Si K) keV. Images are exposure-corrected, 
and darker grey-scales correspond to higher intensities. For the purposes of display, the images 
have been binned into ~5" pixels, and then adaptively smoothed to achieve a minimum S/N = 7. 
J1601 is marked with a cross at the center of the SNR. Image contours of the broadband (1-8 keV) 
image are overlaid in each panel. In (b), the hard feature seen in Fig. [3] is marked with an arrow. 
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Fig. 3.— The 2-8 keV to 1-2 keV hardness ratio map of G330. 2+1.0 as obtained by the XMM- 
Newton EPIC M0S1+M0S2. In the soft band (1-2 keV) map, the 1.4-1.6 and 1.7-1.8 keV bands 
are excluded to remove the bright instrumental lines at ii^ ~ 1.5 and 1.74 keV. For the purposes of 
display, each image has been binned into 15" x 15" pixels, and is adaptively smoothed to achieve 
a minimum S/N = 4. Green image contours are the 1-8 keV image of the SNR as shown in Fig. [2j 
Red contours are the 843 MHz radio image taken from the MOST Supernova Remnant Catalog 
(jWhiteoak &: Green 19961 ) . The angular resolution of the radio image is 43". The position of the 
ceo J 1601 is marked with a cross. The hard feature in the SW of the CCO is marked with a 
white arrow. 
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Fig. 4.— The broadband (1-8 keV) grey-scale image of G330. 2+1.0 obtained by XMM-Newton 
EPIC M0S1+M0S2. The image has been processed in the same way as those in Fig. [2j SW, NE, 
and E regions of the SNR shell are marked with solid lines. Background regions are marked with 
dashed circles. 
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Fig. 5. — The X-ray spectrum from the sheh of G330. 2+1.0. (a) The SW, and (b) the NE shell. 
The best-fit power law model for each regional spectrum is overlaid. The lower panels are the 
residuals from the best-fit model. 
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Fig. 6. — The X-ray spectrum from region E of 0330.2+1.0. The best-fit two component model 
(thick sohd hnes) is overlaid. The dotted and the thin solid lines are the plane-shock (kT = 0.7 
keV) and the power law (T = 2.3) components, respectively. The lower panel is the residuals from 
the best-fit model. 
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Fig. 7. — Radial intensity profiles of the SW shell of G330. 2+1.0 obtained by the Chandra data, (a) 
The bright northern parts (SWl) and (b) the faint southern parts (SW2) of the SW shell, (e) The 
1-7 keV band Chandra ACIS-I image of G330. 2+1.0. SWl and SW2 regions (180" x 40" for each 
region) are marked. The image has been binned into 4" size pixels for the purposes of display. In 
(a) and (b), each regional image has been binned into 8x8 pixels (~4" size), and then is averaged 
over 40" column along the shell to produce a projected 1-D radial intensity profile. The best-fit 
Gaussian model (with a constant underlying background) is overlaid. In (a), the small intensity 
bump just inside of the shell (~25"-45" toward the SNR center) is excluded in the fit. 



